Collective Motion of Moshers at Heavy Metal Concerts 
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Human collective behavior can vary from calm to panicked depending on social context. Using 
videos publicly available online, we study the highly energized collective motion of attendees at heavy 
metal concerts. We find these extreme social gatherings generate similarly extreme behaviors: a 
disordered gas-like state called a mosh pit and an ordered vortex-like state called a circle pit. Both 
phenomena are reproduced in flocking simulations demonstrating that human collective behavior is 
consistent with the predictions of simplified models. 

PACS numbers: 89.75.Kd, 87.15.Zg, 89.65.-s, 47.27.-i 



Human collective behaviors vary considerably with so- 
cial context. For example, lane formation in pedestrian 
traffic [l[, jamming during escape panic and Mexi- 
can waves at sporting events Q are emergent phenomena 
that have been observed in specific social settings. Here, 
we study large crowds (10 2 — 10 5 attendees) of people 
under the extreme conditions typically found at heavy 
metal concerts. Often resulting in injuries [4(, the collec- 
tive mood is influenced by the combination of loud, fast 
music (130 dB [5], 350 beats per minute), synchronized 
with bright, flashing lights, and frequent intoxication Q. 
This variety and magnitude of stimuli are atypical of 
more moderate settings, and contribute to the collective 
behaviors studied here (Fig. [HA)). 

Videos filmed by attendees at heavy metal concerts 
highlight a collective phenomenon consisting of 10 1 — 10 2 
participants commonly referred to as a mosh pit. In mosh 
pits, the participants (moshers) move randomly, colliding 
with one another in an undirected fashion (Fig. [ljB)). 
Qualitatively, this phenomenon resembles the kinetics of 
gaseous particles, even though moshers are self-propelled 
agents that experience dissipative collisions. To explore 
this analogy quantitatively, we obtained video footage, 
corrected for perspective distortions Q as well as cam- 
era instability, and used PIV analysis [H to measure 
the two-dimensional (2D) velocity field on an interpo- 
lated grid. From this data, we calculated the velocity- 
velocity correlation function c vv and noted an absence 
of the spatial oscillations typically found in liquid-like 
systems. Additionally, c vv was well fit by a pure expo- 
nential (R 2 = 0.97) with a decay length of 0.78 m. Taken 
together, these findings offer strong support for the anal- 
ogy between mosh pits and gases. As a further check, we 
examined the 2D speed distribution; previous observa- 
tions of human pedestrian traffic and escape panic led us 
to expect a broad distribution not well described by sim- 
ple analytic expressions (3. fToj| . However, the measured 
speed distribution in mosh pits was well fit by the equilib- 
rium speed distribution of classical 2D gasses (Fig.QJC)), 
otherwise known as the Maxwell-Boltzmann distribution 
[ill ] . These observations present an interesting question: 
Why does an inherently non-equilibrium system exhibit 




FIG. 1. (A) Example image of crowd at a heavy metal con- 
cert. (B) Single video frame illustrating a characteristic mosh 
pit with overlaid velocity field. (C) Speed probability distri- 
bution function (PDF) for the movie in (B) (circles), the best 
fit to the 2D Maxwell-Boltzmann speed distribution (solid), 
and simulated speed distribution (squares). (D) Simulated 
phase diagram plotting the MASHer RMS angular momen- 
tum, demonstrating the existence of mosh pits (gas) and circle 
pits (vortices). (E) Single video frame illustrating a charac- 
teristic circle pit with overlaid velocity field. 



equilibrium characteristics? 

Studies of collective motion in living and complex sys- 
tems have found notable success within the framework of 

Thus, we use a Vicsek-likc 



flocking simulations 
model 



17] to simplify the complex behavioral dynam- 



ics of each human mosher to that of a simple soft-bodied 
particle we dub a Mobile Active Simulated Humanoid, or 
MASHer (SI). Our model includes two species of MASH- 
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ers to reflect the typical crowd at heavy metal concerts 
(Fig. [HA)) Active MASHers are self-propelled, ex- 
perience flocking interactions, and are subject to ran- 
dom fluctuations in the forces they experience. Passive 
MASHers prefer to remain stationary, and are not sub- 
ject to flocking interactions or random forces. Though 
initially both populations are uniformly mixed, we found 
that with sufficient time the flocking-interaction leads 
to a spontaneous phase separation with a dense popu- 
lation of active MASHers confined by passive MASH- 
ers. When noise dominates the flocking and propulsion 
terms, the motion of active MASHers resembles the mo- 
tion in actual mosh pits and quantitatively reproduces 
the experimentally observed statistics (Fig.QJC)). These 
results demonstrate how a non-equilibrium system can 
have equilibrium characteristics: random motions over a 
sufficient time reproduce the statistics of classical gasses 
via the Central Limit Theorem. 

Conversely, when the flocking term dominates active 
MASHcr motion, our model predicts a highly ordered 
vortex-like state [l9| where MASHers again phase sepa- 
rate, but the confined active MASHers move with a large 
non-zero angular momentum (Fig. HJD)). Remarkably, 
this spontaneous phase separation and vortex formation 
is also observed at heavy metal concerts where they are 
conventionally called circle pits (Fig. ffl^E)) Q- In sim- 
ulations, we found an even distribution between clock- 
wise (CW) and counter-clockwise (CCW) motion when 
viewed from above, whereas our observations from con- 
certs show 5% flow CW with the remaining 95% flowing 
CCW (p < 0.001). This asymmetry is independent of 
geographical location, as video data was collected from a 
variety of countries including the United State of Amer- 
ica, the United Kingdom, and Australia. Though the 
origin of this effect is unknown, we speculate it may be 
related to the dominant handedness/footedness found in 
humans. 

The collective behavior described here has not been 
predicted on the basis of staged experiments with humans 
|2Ct [21 [, making heavy metal concerts a unique model 



system for reliably, consistently, and ethically studying 
human collective motion. Further studies in this unique 
environment may enhance our understanding of collective 
motion in riots, protests, and panicked crowds, leading to 
new architectural safety design principles that limit the 
risk of injury at extreme social gatherings. Moreover, 
these concerts have the further advantage of exhibiting a 
rich variety of states such as the wall of death (mosh- 
ers split into two groups separated by an open space 
and, when signaled, simultaneously run at the opposing 
group), collective jumping (locally correlated, globally 
decorrelated), and soft modes [22j in jammed attendees. 
Thus, the extreme conditions at heavy metal concerts 
offer new opportunities for studying a wide range of col- 
lective behaviors that arise in extreme social gatherings 
of large human groups. 



Supplemental Information 

Modeling Moshers with MASHers To model the 
crowd at heavy metal concerts, we simplify the com- 
plex behavioral dynamics of each human mosher to 
that of a simple particle we dub a Mobile Active Sim- 
ulated Humanoid, or MASHer. Analogous to previous 
work 17, [23j], we include the following ingredients in 
our model: (1) Hertzian soft body repulsion (2) self- 
propulsion (3) flocking interactions and (4) Gaussian ran- 
dom noise. The forces on MASHers i are: 
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Here, 



is the center-to-center distance between 



MASHer i and j, ro is the MASHer radius, vq is its 
preferred speed, and Vi is its instantaneous speed. The 
coefficients e, /i, and a are free parameters in our model 
reflecting the relative strength of each term, and the carat 
( * ) denotes unit vectors. The flocking force is a summa- 
tion of the velocities Vj of all the iVj MASHers within a 
radius ra oc k of MASHer i. fj% is a vector whose compo- 
nent are drawn from a Gaussian distribution with zero 
mean and whose standard deviation was varied in the 
simulation. In total we have 7 independent parameters, 
of which all but two are held constant. 

At heavy metal concerts, there are people who do and 
do not participate in moshing. We reflect this in our 
model by having two species. For active MASHers, we 
have, in the appropriate units, vo — I, a varies over the 
interval [0,5], and the standard deviation of the com- 
ponents of ffi varies over the interval [0,3]. For passive 
MASHers Vq = a = \ffi\ = 0. By varying the magnitude 
of the flocking coefficient and the strength of the noise 
for active MASHers, we investigate the phase diagram of 
the model. The remaining parameters are identical be- 
tween the two populations, e = 100, fi = 0.05, ro = 1, and 
mock = 4. 

Simulations We simulate concerts with N = 500 at- 
tendees, of which N a = 150 are active MASHers and 
remaining are passive MASHers. Initially, the active 
MASHers are situated in a circle centered in the mid- 
dle of the simulation box. This choice of initial positions 
is a choice of convenience: any non-zero flocking force 
causes MASHers to self-segregate when started with ran- 
dom initial positions, though this tends to have a large 
timescale (~ 10 3 x (r /v ) time steps) when rflock/ro is 
0(1). Once formed, simulated mosh and circle pits typi- 
cally last for times greater than 10 5 x (ro/uo) time steps, 
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and can switch directions at random intervals 24]. The 
simulation box is of length L = 1.03^/ nr^N » 40.8 to 
give a packing fraction p = 0.94, and is periodic on both 
sides to avoid boundary effects. We use the Newton- 
Stomer-Verlet integration algorithm and cell based neigh- 
bor lists to expedite computation. The phase diagram 
shown in the main text is comprised of 2 x 10 6 individ- 
ual simulations sampled on 10, 000 grid points. For each 
run, we measured the angular momentum about the cen- 
ter of (L/27r)arctan(Im(A)/Re(A)), where 
A = X)t°o 1 ex P ( — ^Trixi/L), with a similar expression 
for Hem- 

Our flocking simulation has three time scales defined 
by the model parameters: the self-propulsion time scale 
Tprop = the flocking time scale Tfl oc k = Vo/a, and the 
collision time scale t co \\ — l/(2rQVop). Both collisions 
and noise tend to randomize motion, whereas flocking 
and self-propulsion tend to homogenize motion. Thus, 
when the time between random collisions t co ii is much 
less than r pr0 p and Tfl oc k, the statistical motion of the 
system is described by a single effective temperature and 
we recover the Maxwell-Boltzmann speed distribution re- 
ported in the main text. 

If we increase the self-propulsion coefficient p (decrease 
''prop)) we find the motion, though random, is no longer 
fit by a Maxwell-Boltzmann distribution. Instead, colli- 
sions between active and passive MASHers on the bound- 
ary of the simulated mosh pit removes energy faster than 
collisions among active MASHers can rethermalize the 
system. Consequently, measurements in silico show a 
radial temperature gradient is established with a higher 
effective temperature at the core of the simulated mosh 
pit and a lower effective temperature at the edge (Fig. [2]). 
The speed distribution in the limit of small r prop is there- 
fore a sum of Maxwell-Boltzmann distributions charac- 
terized by different temperatures. 

Source Code Source code and a phase di- 
agram generating Python script are available 
under the M.I.T. license on github.com at 



https://github.com/mattbierbaum/moshpits. An in- 
teractive Javascript version of the simulation is available 
at |http? //mattbierbaum . github . com/mosh pits . js[ 

The photo in Fig. QJA) was taken and graciously pro- 
vided by Ulrike Biets. JLS and MB also thank D. Porter, 
L. Ristroph, J. Freund, J. Mergo, A. Holmes, A. Alemi, 
M. Flashman, K. Prabhakara, J. Wang, R. Lovelace, P. 
McEucn and S. Strogatz. Field work was independently 
funded by JLS. 




FIG. 2. As the self-propulsion coefficient p increases, we 
find the temperature T(r) has a radial dependence in sim- 
ulated moshpits. This arises from collisions between active 
and passive MASHers that remove energy from the boundary 
faster than collisions among active MASHers can rethermal- 
ize the system. Each line corresponds to values of p starting 
from 0.05 and increasing to 0.50 in steps of 0.05. Simula- 
tions described in the main text correspond to p = 0.05 and 
have essentially a single effective temperature independent of 
radius. 
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